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ABSTRACT 

We present a kinematic diagnostic for peanut-shaped bulges in nearly face-on 
galaxies. The face-on view provides a novel perspective on peanuts which would 
allow study of their relation to bars and disks in greater detail than hitherto pos- 
sible. The diagnostic is based on the fact that peanut shapes are associated with 
a flat density distribution in the vertical direction. We show that the kinematic 
signature corresponding to such a distribution is a minimum in the fourth-order 
Gauss- Hermite moment S4. We demonstrate our method on iV-body simulations 
of varying peanut strength, showing that strong peanuts can be recognized to 
inclinations i ~ 30°, regardless of the strength of the bar. We also consider com- 
pound systems in which a bulge is present in the initial conditions as may happen 
if bulges form at high redshift through mergers. We show that in this case, be- 
cause the vertical structure of the bulge is not derived from that of the disk, that 
the signature of a peanut in s 4 is weakened. Thus the same kinematic signature 
of peanuts can be used to explore bulge formation mechanisms. The observa- 
tional requirements for identifying peanuts with this method are challenging, but 
feasible. 

Subject headings: galaxies: bulges - galaxies: evolution - galaxies: formation - galax- 
ies: kinematics and dynamics - galaxies: photometry - galaxies: spiral 
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1. Introduction 

About 25% of the stellar luminosity in 
the universe comes from the bulges of disk 
galaxies (Persic & Salucci 1992; Fukugita 
et al. 1998). Therefore understanding how 
bulges form is a necessary step in under- 
standing galaxy formation in general. 

If bulges are distinct entities, rather 
than just disk light in excess of an expo- 
nential (van den Bosch et al. 2002; Boker et 
al. 2003), a mechanism for generating them 
separate from disk formation must be con- 
sidered. Because bulges sit at the bottom 
of the potential wells of galaxies, many 
paths for their formation are possible. 
Bulge formation scenarios can be classified 
loosely based on whether the driving mech- 
anism is internal or external. A widely dis- 
cussed example of externally-driven bulge 
formation is in the merger at early times of 
dwarf-sized galactic subunits around which 
disks subsequently grow (Kauffmann et al. 
1993). Observational evidence supporting 
this scenario includes the relatively homo- 
geneous bulge stellar populations in the 
Milky Way (Ferreras et al. 2003; Zoccali et 
al. 2003) and Andromeda (Stephens et al. 
2003), and counter-rotation found in some 
galaxies (Pizzella et al. 2004). 

Discussion of internally driven bulge for- 
mation has focused on the secular evolu- 
tion of disk instabilities. Observational 
evidence supporting secular bulge forma- 
tion includes disk-like, almost-exponential 
light profiles (Andredakis & Sanders 1994; 
Courteau et al. 1996; de Jong 1996; Car- 
current address: Astronomy Department, 
University of Washington, Box 351580, Seattle 
WA, 98195 
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olio et al. 2001; Carollo 1999; Carollo et al. 
1998; MacArthur et al. 2003), occasionally 
disk-like, cold kinematics (Kormendy 1993; 
Kormendy et al. 2002), the correlation be- 
tween the scale-lengths of bulges and disks 
(de Jong 1996; MacArthur et al. 2003) and 
the similar average colors of bulges and in- 
ner disks (Terndrup et al. 1994; Peletier 
& Balcells 1996; Courteau et al. 1996). 
The recent review of Kormendy & Kenni- 
cutt (2004) summarizes our current under- 
standing of, and evidence for, secular for- 
mation of some bulges. 

The bulges of many edge-on galaxies are 
box- or peanut- (B/P) shaped (Burbidge 
& Burbidge 1959; Jarvis 1986). Binney 
& Petrou (1985) constructed axisymmet- 
ric models of B/P bulge systems includ- 
ing cylindrical rotation as observed (Kor- 
mendy & Illingworth 1982). They specu- 
lated that accretion is responsible for cre- 
ating such systems; however observations 
found little evidence of accretion onto them 
(Shaw 1987; Whitmore & Bell 1988, but 
see also Liitticke et al. 2004). A differ- 
ent scenario emerged from 3-D A-body 
simulations, namely formation via secu- 
lar evolution of bars (Combes & Sanders 
1981), either through resonant scattering 
or through bending instabilities (Pfenniger 
1984; Combes et al. 1990; Pfenniger & 
Friedli 1991; Raha et al. 1991). The or- 
bits supporting peanuts have been stud- 
ied extensively (Combes et al. 1990; Pfen- 
niger 1984; Pfenniger 1985; Patsis et al. 
2002a) and shown to generally arise from 
vertically unstable xl orbits. Patsis et 
al. (2002b) showed that these orbits are 
present and peanuts are possible even if 
the non-axisymmetric perturbation is very 
weak. 



2 



Thereafter, observational efforts sought 
to establish the connection between B/P- 
shaped bulges and bars by seeking evidence 
for a bar in edge-on B/P-bulged systems. 
In the case of NGC 4442, the B/P bulge 
is already apparent at an inclination of 
72°, at which the bar also can be recog- 
nized (Bettoni & Galletta 1994). A sec- 
ond such case is NGC 7582 at an inclina- 
tion of 65° (Quillen et al. 1997). In several 
B/P bulges, photometric features of a bar 
have been claimed (e.g. de Carvalho & da 
Costa 1987) but the bar interpretation is 
not unique when only photometry is avail- 
able. The fraction of edge-on bulges having 
B/P shapes is ~ 45% (Liitticke et al. 2000), 
which is consistent with the fraction of 
galaxies containing bars (~ 70%, Knapen 
et al. 2000; Eskridge et al. 2000) once the 
arbitrary orientations of bars to the line- 
of-sight (LOS) are considered. However, 
the most important evidence for the pres- 
ence of bars in B/P bulges comes from a 
comparison of the edge-on gas and stel- 
lar LOS velocity distributions (LOSVDs) 
of iV-body bars (Bureau & Athanassoula 
1999; Athanassoula & Bureau 1999; Bu- 
reau & Athanassoula 2004) and real galax- 
ies (Kuijken & Merrifield 1995; Merrifield 
& Kuijken 1999; Bureau & Freeman 1999; 
Chung & Bureau 2004). 

These edge-on studies have established 
the connection between B/P-shaped bulges 
and bars. However, the degeneracy inher- 
ent in deprojecting edge-on galaxies makes 
it difficult to study other properties of the 
host galaxy. Moreover, while B/P shapes 
are produced by bars, this does not ex- 
clude the possibility that bulges are shaped 
by secular processes, not formed by them. 
Addressing this issue requires an attempt 



at a cleaner separation of bulges, bars and 
peanuts. In face-on systems the viewing 
geometry is well-constrained and bars are 
readily apparent. If we can also recognize 
peanuts in them then we obtain an im- 
portant new perspective on the relation of 
peanuts and bars. For example, this per- 
mits study of the relative sizes of bars and 
peanuts: meager observational evidence 
suggests that these need not be equal (Kor- 
mendy and Kennicutt 2004), in agreement 
with simulations (below). Moreover, for 
inclinations ~ 30°, it becomes possible to 
measure accurately the pattern speed of 
the bar (Debattista 2003) and therefore 
to test for resonances and compare with 
theoretical predictions. It would also al- 
low determination of the fraction of barred 
galaxies with peanuts, which may be dif- 
ferent from the fraction of peanuts with 
bars. And finally, as we will show below, 
the ability to detect peanuts face-on opens 
the possibility of exploring bulge formation 
mechanisms. 

In this paper we examine the kinematic 
signature of peanuts in face-on galaxies. In 
Section 2 we first explore some simple an- 
alytic models to help understand the be- 
havior of more realistic systems. Work- 
ing with Gauss-Hermite moments (Ger- 
hard 1993; van der Marel & Franx 1993), 
we show that the fourth-order LOSVD mo- 
ment, S4, is monotonically increasing with 
g?4, the fourth-order vertical density mo- 
ment. Thus s 4 can be used to probe the 
vertical structure of a disk. We describe 
the iV-body building in Section 3 and in 
Section 4 we present the iV-body models, 
with and without strong peanuts, used in 
this paper and examine their vertical den- 
sity distributions. We show that the main 
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signature of a peanut is in <i 4 , rather than 
in the disk scale-height. We explore the 
vertical LOSVDs of these iV-body models 
in Section 5, showing that s 4 can be used as 
a robust kinematic signature of a peanut, 
independent of bar strength. In Section 6 
we show that moderate inclinations do not 
substantially degrade the diagnostic. Sec- 
tion 7 discusses the required signal-to-noise 
and spectral resolution and our conclusions 
are presented in Section 8. 

2. Exact results 

We first consider some exact models use- 
ful for interpreting the results of iV-body 
simulations. Peanuts constitute a density 
distribution more vertically extended than 
the surrounding disk. Let us denote the 
root-mean-square (RMS) height and verti- 
cal velocity as h z and a z respectively. Devi- 
ations from Gaussian distributions can be 
parameterized by the moments of an ex- 
pansion in Gauss-Hermite functions (Ger- 
hard 1993; van der Marel & Franx 1993). 
The second order term in such an expan- 
sion is related to the RMS. The third order 
term measures deviations which are asym- 
metric about the mean and are therefore 
likely to be small for the vertical density 
and velocity distributions of disk galaxies. 
The fourth order term measures the lowest 
order symmetric deviation from a Gaus- 
sian; it is negative when a distribution is 
broader than Gaussian and positive when 
it is more peaked. We denote the fourth or- 
der Gauss-Hermite moment of the vertical 
density distribution as rf 4 and that for the 
LOSVD as s 4 . Following Gerhard (1993), 
for a vertical LOSVD l(v z ) normalized to 



the projected surface density E, we define 



s 4 = 



l{w)H A {w)e—* w dw (1) 



where w — (v z — vt)jo z and H±(w) = 
-^=(16u> 6 — 48u> 2 + 12). A similar expres- 
sion holds for d 4 . For a particle model, the 
integral becomes a sum and £ is replaced 
by N p , the number of particles in a bin. 

The vertical density extension associ- 
ated with the peanut will correspond to an 
increase in h z (the "scale- height" ) and/or 
a decrease in d 4 . What is the observ- 
able effect on the LOSVD of such varia- 
tions? First consider how o z varies as a 
function of radius in the case where h z is 
constant. In a single-component axisym- 
metric system, the one-dimensional verti- 
cal Boltzmann+Poisson equation is 



d_ 

dz 



ld_ 

pdz 



pv 2 z 



= -AirGp 



(2) 



(e.g. Binney & Tremaine Eqn. 4-38). If 
the system is isothermal, then v z is inde- 
pendent of z (and is therefore equal to a 2 ). 
The solution of Eqn. 2 is 



p(z) = p sech 2 (z/z ) 



(3) 



(Spitzer 1942), where po is the density in 
the mid-plane and 



(J, 



(4) 



(Note that for the isothermal disk, h z = 
z n/^/l2.) Thus if Zq is (nearly) constant 
(as suggested by observations [van der 
Kruit & Searle 1981; de Grijs & Peletier 
1997]) then o\ oc S. Since typically 
H(R) oc e~ R l R ^ then a z oc e - R ^ Rd and it 
is likely that any signature of the peanut in 
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<j z will be swamped by this density-driven 
radial variation. 

Next we use simple models to explore 
the signature in the LOSVD of a flat- 
topped vertical density distribution. Con- 
sider Camm's (1950) series of analytic so- 
lutions of the collisionless Boltzmann equa- 
tion for systems stratified in plane-parallel 
layers of infinite extent. This is a good lo- 
cal approximation to real galaxies at low 
z when the rotation curve is flat (van der 
Kruit & Freeman 1986). In Camm's model 
III, the density distribution is given by 
p(z) = pocos 2 ~ 2 / n 9, where n > §, and the 
parameter 9 is defined by the relation 

f e 

z = A sec 1 " 27 " d<P (5) 
Jo 

with A some constant. Several examples of 
this density distribution for different n are 
presented in Figure 1. The corresponding 
distribution function is 

f(w,$) = C{n){2ncos 2/n 9 -w 2 ) n ~l, (6) 

where C{n) is a normalization constant. 
These densities and distribution functions 
can be integrated numerically to compute 
d,4 and S4. The results are presented in 
Figure 2. The limit n — > 00 corresponds to 
the isothermal sheet, in which the distri- 
bution function has the same Gaussian de- 
pendence on velocity at all heights. Thus 
in this limit, S4 = 0. On the other hand, 
the isothermal sheet has a sech 2 z profile 
which is more peaked than a Gaussian and 
therefore has d± > 0. At smaller n, the 
density profile becomes increasingly flat- 
topped leading to rf 4 < which drives 
S4 < 0. Two properties of S4 make it an 
excellent probe of <i 4 . First is the fact that 
s 4 increases monotonically with d 4 , which 



makes S4 an observable surrogate for the 
unobservable d 4 . Second, s 4 ^ d 4 , so that 
the vertical velocity distribution is gener- 
ally broader than the density distribution, 
which makes it observationally robust. 
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Fig. 1.— Model III of Camm (1950). Vari- 
ous vertical density profiles, with n increas- 
ing in order of increasing maximum z, from 
n — 1.6 in steps of 1, are shown. 

3. iV-body Systems 

We use iV-body models with different 
initial conditions and evolved on different 
codes to study the kinematic signatures 
of face-on peanuts. Since gas is gener- 
ally depleted within bars, and moreover 
dissipates its vertical energy, it is not a 
good tracer of face-on peanuts. We there- 
fore focus on only the stellar kinematics of 
peanuts, and our simulations are all colli- 
sionless. Table 1 lists all the simulations 
used in this study. 

3.1. Rigid-halo simulations 

The highest resolution simulations in 
this paper were run on a 3-D cylindrical 
particle-mesh (PM) grid code (described 
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n d 4 

Fig. 2. — The Gauss- Hermite moments d 4 
and S4 of Camm's (1950) model III. The 
left panels show the variation of these pa- 
rameters with n, while the panel on the 
right plots s 4 versus d 4 . The dashed diag- 
onal line indicates s 4 = d± 

in Sellwood & Valluri 1997). The main 
advantages of this code for the present 
study are that it permits high spatial and 
mass resolutions; indeed these simulations 
used > AM disk particles and force soften- 
ing e < 2^d/3, where z d is the (constant) 
Gaussian vertical scale-height of the initial 
conditions. Since Gauss-Hermite moments 
generally require high signal-to-noise ratios 
to be measured reliably, a large number of 
particles is desirable. At the same time, 
the high force resolution ensures that the 
vertical motions of particles are well re- 
solved. 

The rigid halos were represented by ei- 
ther a spherical logarithmic potential, 

$ L (r) = | Hr 2 + rl), (7) 



or a Hernquist (1990) model 



where r h is a halo scale-radius, t>h is a char- 
acteristic halo velocity and M h is a halo 
mass. We define t>h = y/GM h /r h for the 
Hernquist halos. 

The initially axisymmetric disks were all 
Sersic type 

Pd (R,z) oc (1 - /bjMe-W^'V*^ 2 

(9) 

where /b is the fraction of the active (i.e. 
bulge+disk) mass which is in the bulge, M 
is the active mass, R d is the disk scale- 
length, z d is the Gaussian thickness and 
n is the Sersic index (n — 1 correspond- 
ing to an exponential profile and n = 4 
to a de Vaucouleurs profile). Disk kine- 
matic setup used the epicyclic approxima- 
tion with constant Toomre-Q and the ver- 
tical Jeans equation to set vertical motions 
appropriate for a constant Gaussian thick- 
ness. We use units where R d = M = G = 
1, which gives a unit of time (R^/GM) 1 ^ 2 . 

Bulges were generated using the method 
of Prendergast & Tomer (1970), where a 
distribution function is integrated itera- 
tively in the global potential until conver- 
gence. We used the isotropic distribution 
function of a lowered polytrope, truncated 
at r b 

/(f,t7)oc[-E(f,t;)] 1/2 -[-E max ] 1/2 . 

(10) 

Here E max = $ to t(r b ), the total potential 
at r b in the disc plane. For all bulges we 
set rb = 0.78. The bulges in runs B2 and 
B3 were fully rotating, while that in run 
Bl had no rotation. Further details of the 
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Strong 
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Table 1: The sample of simulations used in this paper. TV* is the number of disk+bulge 
particles, n is the index of the initial Sersic disk, and e are the scale-length, scale- 

height of the initial disk and softening length, Q is the initial disk Toomre-Q and is 
the bulge mass as a fraction of the total (disk+bulge). In column "Halo" we describe the 
type of halo used: logarithmic, Hernquist or NFW. r h and v h are the halo scale-length and 
characteristic velocity, respectively. In column "Peanut" we give a qualitative description of 
the peanut: strong, weak or none. For the live-halo system LI, we give the minimum Q, 
is for a sech 2 profile, r h = r vir and = V vir in kpc and km s -1 respectively. 



compound system setup method used can 
be found in Debattista & Sellwood (2000). 

The polar grids were N R x x N z = 
60 x 64 x 225 or larger. For all the PM sim- 
ulations, the vertical spacing of the grid 
planes, 5z, was set to 0.0125. We used 
Fourier terms up to m = 8 in the po- 
tential, which was softened with the stan- 
dard Plummer kernel. Time integration 
was performed with a leapfrog integrator 
with a fixed time-step, 5t = 0.02 for simu- 
lations B1-B3, 5t = 0.0025 for run R4 and 
St = 0.01 for all the rest. 

3.2. Live-halo simulation 

The disadvantage of the PM code is that 
we needed to use a rigid halo. There- 
fore in run LI we used a lower mass res- 



olution live-halo simulation run with pkd- 
GRAV (Stadel 2001), a multi-stepping, par- 
allel treecode. 

The live-halo model was built using the 
technique developed by Hernquist (1993; 
see also Springel & White 1999). We start 
with an isotropic NFW halo (Navarro, 
Frenk & White 1996) with virial radius, 
R vir = 108 kpc, circular velocity at the 
virial radius, V vir = 76 km s _1 and virial 
mass M vir = 1.5 x lO n M . Then adia- 
batic contraction of the halo due to the 
presence of the disk is taken into account 
assuming that the spherical symmetry of 
the halo is retained and that the angular 
momentum of individual dark matter or- 
bits is conserved (see Springel & White 
1999). The disk mass fraction relative to 
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the halo virial mass, f d = M d /M vir = 0.08. 
We used an exponential disk with scale- 
length Rd = 1.99 kpc and a sech 2 (z/zd) 
vertical profile. We set Zd = 0.15i?d and 
softening length e = 50 pc. The veloc- 
ity field of the disk was calculated as in 
Springel & White (1999) assuming the ra- 
dial and vertical velocity dispersions are 
equal, <jr = &z, with <jr chosen to give 
minimum Toomre-Q = 1.2. Then the az- 
imuthal velocity dispersion is determined 
from or using the epicyclic approximation. 

4. Vertical Density of the iV-body 
Models 

The evolution of model Bl has been de- 
scribed in Debattista (2003), while runs B2 
and B3 formed part of the preliminary sur- 
vey for the Milky Way modeling described 
in Bissantz et al. (2004). The evolution of 
most of the remaining models will be de- 
scribed elsewhere (Debattista et al. 2005 
in progress). Here we are interested pri- 
marily in the final systems not in details of 
their evolution. Except for run R3, which 
formed only a very weak oval distortion, all 
these simulations formed bars. The edge- 
on view of the simulations R1-R8 is pre- 
sented in Figure 3. Throughout this pa- 
per, we use a convention where the bar 
is along the rr-axis and the z-axis is per- 
pendicular to the disk. Runs R1-R4 con- 
tain prominent peanuts which were pro- 
duced by bending instabilities (Raha et al. 
1991; Merritt & Sellwood 1994). To bet- 
ter present these peanuts, in Figure 4 we 
present the edge-on projected density of 
particles in the narrow range —0.5 < y < 
0.5, i. e. we show only a narrow slice of each 
model extending to about the minor-axis 
of the bar. This gives a better apprecia- 



tion of the peanuts which will be sought 
in the face-on view, where the disks do 
not mask peanuts. Peanuts can form in 
weakly barred systems (Patsis et al. 2002); 
the peanut in run R3 formed in the pres- 
ence of only a very weak oval. As a re- 
sult, this peanut is almost axisymmetric. 
Run R5 contains a weak peanut while runs 
R6-R8 have no peanuts at all. Run R8 
is identical to run Rl except for one im- 
portant detail: we forced symmetry about 
the mid-plane. Therefore, although the 
velocity ellipsoid is very anisotropic with 
a z /a R ~ 0.25 through most of the bar, the 
system could not bend and did not develop 
a peanut. 

Comparing the peanuts visible in Fig- 
ures 3 and 4 with the bars seen in the face- 
on surface density overlaid on Figure 5, it 
is apparent that the peanut in run Rl is 
smaller then the bar. Kormendy and Ken- 
nicutt (2004) noted that the peanut in the 
moderately inclined galaxy NGC 7582 was 
significantly shorter than the bar, and wor- 
ried that this may be a problem for a sec- 
ular peanut formation scenario. Run Rl 
shows that peanuts need not fill the entire 
major axis of a bar. In contrast, the peanut 
in run R2 extends to about the ends of the 
bar. 

Figure 5 shows maps of h z for runs Rl- 
R8. On the bar's major axis h z increases 
radially outwards, reaching a maximum 
when a peanut is present. In most cases, 
the local maximum in h z occurs close to 
the peanut. On the other hand, in runs 
R6-R8, none of which contain a peanut, h z 
increases throughout. 

The maps in Figure 6 show <i 4 for runs 
R1-R8. Minima in <i 4 correlate well with 
the location of peanuts and there are no 
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significant minima in d± other than at R = 
in the absence of a peanut. 

For better comparison of the different 
models, Figure 7 shows the density, h z and 
di profiles along the major axes of the face- 
on bars. 

Run R1 Run R5 

Run R2 Run R6 

Run R3 Run R7 

Run R4- Run R8 




I -2 2 -2 

x x 



Fig. 7. — A summary of the face-on den- 
sity moments along the bars' major axis 
for runs R1-R8. 

4.1. h z versus d 4 

Are the second and fourth density mo- 
ments equivalent ways of defining the 
peanut? Figure 7 shows that the maxima 
in h z corresponds very well to the minima 
in G? 4 in runs Rl and R2. In run R3 the 
peak of h z is at a different radius from 
the minimum of d±, while in run R4 h z 
and c/4 have maxima at about the same 
point. Thus the maximum in h z and the 
minimum in d^ are not equivalent ways of 
defining the peanut. 



Which of h z and d± is the better tracer 
of peanuts? Figure 7 shows that the 
major- axis h z profiles of runs Rl and R5 
are rather similar, suggesting very similar 
peanuts; Figures 3 and 4 show that this is 
far from being the case. In contrast, their 
c?4 profiles in Figure 7 are very different. 
For this reason c? 4 is a better measure of 
the presence and strength of a peanut. The 
reason why h z is not an optimal peanut di- 
agnostic is that it is partly determined by 
the local projected density, as suggested 
by the correlation, evident in Figure 7, be- 
tween the depth of the central minimum 
in h z and the central concentration. 

5. Vertical Kinematics of Peanuts 

5.1. The absence of a peanut signa- 
ture in a z 

Figure 8 shows maps of a z for runs Rl- 
R8. No sign of the peanuts is evident in 
these maps. We found that, in the region 
—2.5 < x, y < 2.5 of all the models, a z cor- 
relates very strongly with S, even off the 
bar's major axis. The signature of a peanut 
is buried in the small scatter in a z at fixed 
S, making the peanut hard to distinguish 
from a z . 

5.2. A peanut signature in a^/S 

Following the discussion in Section 2, 
it is unsurprising that we cannot identify 
peanuts from o z profiles. Since a 2 z ~ £ in 
the isothermal disk (Eqn. 4), it is worth 
exploring whether cj^/E is better than o z 
at locating peanuts. In Figure 9 we plot, 
for runs R1-R8, cr^/S on the bars' major- 
axes. In all cases, the profiles of logt^/E) 
are rather similar to those of h z (although 
there is not a one-to-one correspondence). 
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Fig. 3. — The edge-on view of runs R1-R8. In all cases, the bar has been rotated into the 
x-axis, so it is viewed side-on. 

Fig. 4. — The edge-on view of the systems as in Figure 3 but with \y\ < 0.5 only shown. 

Fig. 5. — The face-on view of runs R1-R8 showing the color-coded RMS height (h z ) of 
particles. The contours show the projected surface density, S. 

Fig. 6. — The face-on view of runs R1-R8 showing the fourth-order Gauss-Hermite moment 
of the vertical density distribution, d 4 . 



Fig. 8. — The face-on view of runs R1-R8 showing o z . 



Run R1 Run R5 

Run R2 Run R6 

Run R3 Run R7 

Run R4 Run R8 




I -2 2 -2 



Fig. 9. — A summary of the face-on kine- 
matic moments along the bar's major axis 
for runs R1-R8. The peanut is not appar- 
ent in the velocity dispersion profile but is 
prominent in the s 4 profile. 

However, a comparison of Figures 7 and 9 



shows that the minimum in <i 4 , which we 
showed above to be a good indicator of a 
peanut, does not correspond to any special 
point in the cr^/E profile. Therefore of/ 1! 
has the same limitations peanut diag- 
nostic as h z , although it may still be use- 
ful in distinguishing between peanut and 
peanutless systems. 

5.3. The peanut signature in s 4 

Figure 10 plots s 4 for runs R1-R8. A 
qualitative difference between the peanut 
systems and the peanutless ones is evi- 
dent - - two negative minima at the lo- 
cation of the peanut on the bar's major 
axis if a peanut is present. In runs Rl- 
R4 there is a considerable variation in bar 
strength, but in each case, the negative s 4 
minimum criterion recognizes the peanut, 
demonstrating that it does not depend on 
bar strength. 

In run R5, which produced a weak 
peanut, s 4 remains greater than zero on 
the bar major-axis and the only minimum 
is at R = 0. Thus the kinematic diagnostic 
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cannot identify weak peanuts. Other than 
at R = 0, no significantly negative min- 
imum in S4 occurs in runs R6-R8 which 
lack a peanut. 

Figure 11 plots S4 versus cLt on the major 
axes of the bars. As predicted by the ana- 
lytic models of Section 2, the minimum in 
g?4, which is a tracer of the peanut, corre- 
sponds to the minimum of S4. Therefore 
minima in S4 are an excellent kinematic 
peanut diagnostic in these face-on systems. 

5.4. Live-halo simulation 

Mainly because of the lower mass reso- 
lution, which results in significantly lower 
S/N in the Gauss-Hermite moments, we 
used the live-halo run LI only to confirm 
results of the rigid-halo simulations. De- 
spite the lower S/N, we were still able to 
identify clear and well-matched minima in 
g?4 and S4 in the region of the peanut. In- 
deed the properties of the peanut in LI are 
very similar to those in run Rl although 
the initial conditions, including vertical 
structure, were quite different. This gives 
us confidence that the peanut diagnostic 
developed from rigid-halo simulations is 
not an artifact of the rigid halos. 

5.5. Simulations with bulges 

Now we consider simulations which in- 
clude a bulge in the initial conditions. In 
Figure 12 we plot the edge-on views of 
these systems. A weak peanut is present in 
run Bl, which is masked in the full edge- 
on view. A stronger peanut is present in 
run B2 and an even stronger one in B3. 
Comparing with Figures 3 and 4, the pure 
disk components of runs Bl, B2 and B3 
are most like runs R5, R4 and Rl respec- 
tively, albeit only approximately. Figure 



13 plots maps of and s 4 and the major 
axis profiles are presented in Figure 14. In 
all three runs, the slope of the profile of 
h z has a break which allows the peanut to 
be recognized. Thus the peanuts are still 
evident in cx^/S, including again the weak 
peanut in Bl. All three models have min- 
ima in their <i 4 profiles, but those in runs 
Bl and B2 are broad ones extending down 
to R — 0. In S4, no minima are visible in 
run Bl (and thus no peanut is identified, 
as was the case also for the weak peanut 
in run R5), while the usual minima iden- 
tifying a peanut are clear in run B3. The 
case of run B2 is more interesting. If the 
bulge in this system were dark, a peanut 
would stand out clearly in s 4 (cyan line); 
with the addition of the bulge, the differ- 
ent bulge kinematics, especially the signif- 
icantly higher a z , perturb the net s 4 , hid- 
ing the presence of a peanut in the sense 
that only a single broad minimum down to 
R = remains in its profile. 

Fig. 12. — The disk+bulge simulations. 
Panels (a) show the bulges, panels (b) 
the disks and panels (c) the bulges+disks. 
Panels (d) and (e) show the bulges+disks 
and disks in the range \y\ < 0.5. 



Fig. 14. — The photometric profiles of runs 
B1-B3 on the bars' major axis. Red (dot- 
ted) is the bulge, cyan (dashed) is the disk 
and black (solid) is disk+bulge. 

6. The Effect of Inclination 

Exactly face-on galaxies are rare; the 
probability for a galaxy to be within 5° 
of edge-on is over 20 times larger than for 
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Fig. 10- 



The face-on view of runs R1-R8 showing S4. 
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Fig. 11. — The variation of S4 versus d± on the bar's major axes for runs R1-R8. 



Fig. 13. — The bulge simulations B1-B3. 



it to be within 5° of face-on (8.7% versus 
0.4%). The inclination needs to be within 
24° of face-on before its probability is equal 
to that for 5° from edge-on. A sample of 
exactly face-on galaxies may therefore be 
hard to obtain. Thus it is necessary to ask 
what happens to the kinematic signature of 
a peanut when it is viewed not quite per- 
fectly face-on. Is it possible that the neg- 
ative minimum in S4 signature of a peanut 
is erased, or induced where no peanut is 
present, for other than an exactly face-on 
orientation? 



Once a system is no longer perfectly 
face-on, the observed LOS velocity dis- 
persion, (7i os , relative to which s 4 is de- 
fined, includes contributions from the ra- 
dial ((Jr) and tangential (cr^) dispersion 
components. If % is the inclination angle, 
is any angle in the disk's plane measured 
relative to the inclination axis, a = cr^/a R 
and (3 = o- z /cr R , then the contribution of 
cr z to <7i os is 



O"los 



(3 2 sin 2 sin 2 i + a 2 f3~ 2 cos 2 sin 2 i + cos 2 i 

(11) 
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Run R1 Run Run B2 Run B3 
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Fig. 15. — A comparison of the density 
profiles along the bar's major axis in runs 
Rl and B1-B3. The central concentration 
is largely due to initial conditions in runs 
B1-B3 but results from secular evolution in 
run Rl. The dashed black line shows the 
scaled-profile along the bar major- axis of 
NGC 4477. 

This has a maximum along the disk's ma- 
jor axis (0 = 0) and a minimum on the 
minor axis (0 = 90°). A crude estimate 
of <7 z /a\ os can be obtained assuming that 
a 2 = 1/2 (i.e. a flat rotation curve). Then 
if (5 — 0.293 (the minimum value required 
for stability [Araki 1985]), we estimate that 
the contamination may be as high as 10% 
already at an inclination of 10°. However, 
(3 = 0.293 is extreme: in the Solar neigh- 
borhood, Dehnen & Binney (1998) find 
P = 0.53 ±0.07 while Gerssen et al. (2000) 
find even larger {3 in earlier Hubble types. 
The contamination of o~ z in <7i os is still less 
than 10% at i = 30° on the disk's major 
axis if j3 = 0.53. 

Since these estimates are based on sim- 
plifying approximations, we also explored 
the effect of inclination directly on the N- 
body simulations to i = 40°. Figure 16, 
presents S4 on the bar major-axis for runs 
R1-R8 inclined at i — 30° and with the bar 



oriented at 0° < < 90°. 

For strong peanuts, the negative S4 min- 
imum criterion still distinguishes between 
peanut and peanutless systems up to an 
inclination of ~ 30°. When a peanut 
is present, the two minima on opposite 
sides of the bar become asymmetric as 
increases to 90°. Therefore ^ 45°, 
is a more favorable orientation for find- 
ing peanuts. In some instances, inclina- 
tion produces negative minima in S4 off 
the bar's major axis in both peanut and 
peanutless systems. However, inclination 
leads to, at most, only shallow minima 
in peanutless systems, although an over- 
all negative s 4 can result when no peanut 
is present. (Thus a negative S4 without 
a minimum is not by itself sufficient as a 
peanut diagnostic.) 

In run R8 we prevented bending by forc- 
ing mid-plane symmetry, which resulted in 
a final <j z /<Jr ~ 0.25. When this system is 
viewed at an inclination of 30°, two shal- 
low minima in S4 appear at all (see Figure 
16) even though no peanut is present. As 
expected from the analytic estimate above, 
inclination has a much larger effect on the 
vertical kinematic moments of vertically 
cold systems. As this is an unrealistically 
anisotropic system (symmetrization about 
the disk plane having inhibited the bend- 
ing instability), this simulation represents 
an extreme extent to which inclination in- 
troduces minima in S4 when no peanut is 
present. 

Two animations accompany the online 
version of this paper. These show the ef- 
fect of inclination on the S4 moments of 
simulations Bl and B3. 
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Fig. 16. — The s 4 profiles on the bar 
major- axis seen face-on (black line) and at 
i = 30° for runs R1-R8. The bar makes 
an angle = 0°, 45° and 90° for the cyan, 
green and red lines, respectively. In our 
notation x > is the side nearer the ob- 
server. 

7. Observational Requirements 

Measuring kinematic Gauss-Hermite 
moments requires high S/N spectra (Ben- 
der et al. 1994). In our iV-body measure- 
ments, this has been possible because of 
the large number of particles. By resam- 
pling experiments, we determined that a 
S/N ^ 50 is required to measure S4 suffi- 
ciently accurately to identify a peanut, in 
good agreement with Bender et al. (1994). 
Fortunately bars are generally bright fea- 
tures which helps improve the S/N. 



The kinematic signature of a peanut 
is strong for some distance on the bar's 
minor-axis (Figure 10). Thus very precise 
placement of the slit along the bar's ma- 
jor axis is not necessary. This also allows 
the widest slit consistent with the neces- 
sary spectral resolution. 

The spectral resolution, R, required de- 
pends on the value of a z which varies from 
galaxy to galaxy. The Milky Way has 
o z ~ 100 km s" 1 (Kuijken 2003). Thus one 
would need R ~ 2500 to find a peanut in 
a face-on galaxy like the Milky Way. If, on 
the other hand, o z ~ 30 km s _1 (e.g. Bot- 
tema 1993), R ~ 8500 would be needed. 

8. Discussion 

B/P-shaped bulges are common in edge- 
on disks and simulations show that these 
can form by secular evolution in barred 
galaxies. In this paper we explored the sig- 
nature of peanuts on the vertical density 
and the resulting kinematics. We showed 
that both h z and are affected by a 
peanut, but the two are not equivalent sig- 
natures of a peanut. The preferred peanut 
signifier is g? 4 since it distinguishes between 
weak and strong peanuts. 

The vertical velocity dispersion, a z is a 
poor diagnostic for peanuts because it de- 
pends on the local density. The quantity 
cr^/S, which in an isothermal disk would 
trace h z , factors out some of this depen- 
dence and is able to identify peanuts, even 
weak ones, at breaks in its slope. However 
it is unable to quantify peanut strength 
and correlates poorly with d 4 . This pa- 
rameter may also be prone to systematic 
effects from variations in mass-to-light ra- 
tios if, for example, the bulge and the disk 
are composed of different stellar popula- 
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tions, as would happen if the bulge formed 
at high redshift through mergers. Never- 
theless, this is a useful peanut diagnostic 
that is worth testing in real galaxies. 

An excellent kinematic diagnostic of 
face-on peanuts are negative double min- 
ima in the Gauss- Hermite moment S4. The 
negative s 4 minimum signature of a peanut 
holds for any bar strength down to the 
weakest of ovals and may therefore be 
used to search also for peanuts in unbarred 
galaxies. This diagnostic is not too sensi- 
tive to inclination for % < 30°, with nega- 
tive minima on a bar's major-axis contin- 
uing to be associated with the presence of 
a peanut. However, inclination leads to an 
asymmetry between the two sides of the 
bar as the bar orientation approaches the 
minor axis. Thus a bar oriented within 
~ 45° of the line-of-nodes is ideal for a 
peanut search. 

8.1. Bulge formation mechanisms 
and vertical structure 

Two competing models of bulge forma- 
tion — via internal secular evolution and 
via external drivers — each account for a 
significant body of observational evidence, 
suggesting that both processes play some 
role (e.g. Wyse 2004). Thus we need to ask 
which process dominates in which galaxies. 
The results here suggest a novel observa- 
tional program to address this question by 
targeting the degree of decoupling of the 
vertical structure of bulges and bars. To 
be concrete, consider runs R4 and B2, both 
of which contain a peanut. Since the bulge 
is already present in the initial conditions 
in run B2 (as would be a bulge formed in 
an early merger), the presence of a peanut 
cannot be used to address whether the 



bulge formed by secular evolution or not. 
Photometrically the two have similar den- 
sity profiles (Figure 15) since run Rl ac- 
quired a central density concentration by 
secular evolution (Hohl 1971; Debattista 
et al. 2004). These profiles are typical of 
the major-axis profiles of real galaxies: the 
dashed line in this figure shows the J-band 
profile of the nearly face-on galaxy NGC 
4477, taken from the online 3 Frei catalog 
of galaxies (Frei et al. 1996). However the 
vertical structure of the bulge in run B2 is 
qualitatively different from that in run Rl 
because it is not derived from the disk's. 
Because of this, neither the d 4 nor the s 4 
profiles have the kind of separated minima 
associated with peanuts in which the cen- 
tral concentration forms purely by secu- 
lar evolution. Peanuts have been shown 
to be visible in 45% of edge-on galaxies 
which means that they are even more com- 
mon since projection hides some fraction of 
peanuts. Thus a kinematic survey of face- 
on barred galaxies should turn up a large 
fraction of galaxies with the negative s 4 
minima signature of a peanut if bulges are 
built largely by secular evolution of bars. 
On the other hand, if bulges formed largely 
through mergers of dwarfs then peanuts 
would need to be stronger to be identifi- 
able in face-on kinematics. 

8.2. The effect of gas 

The simulations presented here were all 
collisionless. If gas funneled by a bar plays 
an important role also in bulge formation, 
naively it would seem likely that the ver- 
tical kinematic signature of a peanut be- 
comes confused. However, in a barred po- 

3 Available at http://www.astro.princeton.edu/~frei/catalog.htm 
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tential, gas sinks to small radii where its 
kinematics and of stars formed from it do 
not perturb those from the peanut fur- 
ther out. Lower resolution hydro+iV-body 
live-halo experiments we have run and 
which will be presented elsewhere show 
that peanuts can still be recognized in that 
case with the same stellar kinematic diag- 
nostic. Therefore an observational sur- 
vey of face-on barred galaxies to look for 
peanuts appears worth undertaking. 

We would like to thank the anonymous 
referee for comments which helped improve 
the presentation of this paper. V.P.D. 
thanks Enrico Maria Corsini, Sven De Ri- 
jcke and Ortwin Gerhard for fruitful dis- 
cussion. 
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